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ABSTRACT 

We discuss the problem of how to calculate the distance between two cosmological 
objects given their redshifts and angular separation on the sky. Although of a funda- 
mental nature, this problem and its solution seem to lack a detailed description in the 
literature. We present a new variant of its solution and quantitatively assess the most 
commonly used approximation. 
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1 INTRODUCTION 

In cosmology and extragalactic astronomy one frequently 
needs to calculate the distance between two objects given 
their redshifts and their angular separation on the sky. As 
larger and larger cosmological volumes are probed by wide 
field redshift surveys such as 2dF and the Sloan Digital Sky 
Survey, or by QSO absorption line studies, the effects of 
non-Euclidean geometry become increasingly important. In 
the past, authors have frequently relied on approximations 
when calculating the distance between two cosmological ob- 
jects, presumably for calculative ease. Such approximations 
are valid only for small distances and are particularly use- 
ful when examining or highlighting the feasibility of geo- 
metrical mea ns to measure cosmological p arameters as wa s 
don e e.g. by |Alcock fc Paczyriski (197S| ), [Phillipps (1993 ) 
and Popowski et al. (1995 ). However, many practical appli- 



cations are not limited by computing time and since an ex- 
act solution to the problem exists the approximations seem 
unnecessary. These applications include the construction of 
the real-space two-point cor relation function of various ob- 
jects such as galaxies (e.g. [Yoshii, Peterson, & Takahara 



19931) and QSO absorption systems (Williger et al. 1996^ 



Dinshaw & Impey 1996), as well as studies of the effect of 



local sources of ionizing radiation on their surroundi ng in- 
tergalactic medium (e.g Fernandez-Soto et al. 1995). The 
latter application actually requires knowledge of not only 
the distance between two objects but also of the redshift ex- 
perienced by a photon travelling from one object to another. 

Given the fundamental nature of this problem we feel a 
detailed, explicit treatment is called for. In this short article 
we present a new variant of the solution to the distance 
problem (Sections |^ and |§]), discuss its relation to existing 
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approaches (Section^) and finally investigate the validity of 
the most frequently used approximation (Section 

For clarity and brevity we have limited ourselves in this 
paper to homogeneous Friedmann (zero-pressure) cosmolo- 
gies with no cosmological constant (A = 0). The inclusion 
of A renders some of the explicit expressions non-analytical 
and thus (in the context of this paper) unnecessarily com- 
plicates matters. 



2 THE DISTANCE BETWEEN ANY TWO 
OBJECTS 

We begin by writing the familiar Robertson- Walker line el- 
ement as: 

ds 2 = -cW + a{t) [d X 2 + E 2 (x)(d6» 2 + sin 2 9 d(f> 2 )] , (1) 
where 




k = +1 
k = 
k = -1. 



(2) 



Putting the Earth at the origin of the coordinate system one 
can use this metric and the Friedmann equations to calculate 
the distance from Earth (at x = 0) to an object at redshift 
z, corresponding to a comoving coordinate Xj 

r = a E[x(«)] 
c 1 



H ql 1 + z 

x q z + (qo - 1) ( y/T+ 2q z - 1 



(3) 



(e.g. Misner, Thorne, & Wheeler 1973), where ao, -Ho and 



go are the scale factor, Hubble and deceleration parameters 
at the present epoch (subscript 0). 
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Now consider an object 1 (the 'receiver') observed on 
Earth today at Z\ and an object 2 (the 'emitter') at z 2 sep- 
arated by an angle a on the sky (cf. Fig. Q). Object 2 emits 
a photon towards object 1 which is received by object 1 at 
the same time as object 1 emits the photon we receive from 
it today (i.e. at the epoch corresponding to Zi). What is the 
distance, r' 2 , between these two objects at the time of the 
photon reception and what is the redshift, z' 2) of the photon 
as observed by object 1? 

An observer at object 1 would write equation (^|) as 

r' 2 = OiE(xi) 
c 1 



Earth 



qiz' 2 + (qi - 1) ( y/l + 2qiz' 2 - 1 



(4) 



where x'2 ls the comoving coordinate distance between ob- 
jects 1 and 2 (cf. Fig. [I]) and 



ai = 



Hi 



0(1 



l + Zl 



H (l + zi) v /l + 2q zi, 
1 + 21 



(5) 



go: 



1 + 2q Z! 

are the scale factor, Hubble and deceleration parameters at 
the time object 1 emitted the photons we receive today. 
Thus we see that the problems of calculating r' 2 and z' 2 are 
equivalent since knowledge of one provides knowledge of the 
other via equation (^). Here we choose to find z' 2 . 

Note that r 2 does not describe in general the shortest 
distance between objects 1 and 2 along a t — const hypersur- 
face of spacetime. Nevertheless, in many applications r' 2 is 
the quantity of interest. E.g. when considering the radiative 
effect of a QSO on a nearby object one needs the luminosity 
distance between the two which is given by r' L2 — r' 2 (l + z' 2 ). 
In any case, it is most practical to solve the problem for z 2 
and then calculate the distance required for a given appli- 
cation from z 2 . 

As we are dealing with three points (Earth and objects 
1 and 2) in a 3-dimensional space (described by the spatial 
part of the metric, equation ^) it seems intuitive that it must 
be possible to reduce the problem to two dimensions as one 
can always find a 2d hypersurface that contains all three 
points. Since we are interested in measuring 3d distances 
the hypersurface should be chosen such that the distance 
between any two points of the surface (as measured along 
the surface) is identical to the 3d distance between the same 
two points. We call such surfaces totally geodesic. It is clear 
that a totally geodesic hypersurface containing a given set 
of three points can only be constructed from the geodesies 
connecting the three points. Since these geodesies are unique 
there can be only one such surface. 

In equation (|l|) we introduced a polar coordinate system 
(x,0, </>)• Clearly, the hypersurface described by <f) — const 
is totally geodesic. (Note that the surfaces \ — const and 
9 = const are not.) Since the curvature of the 3d space un- 
der consideration is constant, one can generate all totally 
geodesic hypersurfaces from any given one by mere trans- 
lations and rotations. Therefore, for a given set of three 
points there must exist a coordinate system (x, 8, <t>) such 
that (j> = const describes the unique totally geodesic hyper- 
surface containing these three points. Since this new coordi- 




Figure 1. In the case k = +1 objects 1 and 2 and Earth form 
a triangle with geodesic sides (thick lines) on the surface of a 2d 
sphere of radius a. The Earth is at the origin of the coordinate 
system ('north pole'). From the centre of the sphere the angle 
between Earth and object 1 is jq. A photon emitted by object 
1 towards Earth travels along the geodesic connecting the two 
(thick line). The length of this path (= distance between object 1 
and Earth) is a\i- aE(xi) = asinxi is the distance from object 
1 to the central axis of the sphere. Essentially this is the angular 
diameter distance or luminosity distance (modulo factors of 1 + z) 
from Earth to object 1. \' 2 is the unknown angle between objects 
1 and 2 at the centre of the sphere and a\' 2 is the unknown length 
of the connecting geodesic. 



nate system can be constructed from the old one by transla- 
tion and rotation, the form of the metric in this new system 
is identical to equation ([!]). Restraining this metric to the 
<f> — const hypersurface we have 

dl 2 = a 2 [d X 2 + E 2 ( X )d<9 2 ] . (6) 

Thus we can see that the triangle Earth-object 1-object 
2 lies either on a 2d sphere, a plane or a 2d hyperboloid 
(k = +1,0, —1) embedded in 3d Euclidean space. This tri- 
angle has geodesic sides XiiX2 5 and X2 (connecting the two 
objects) and the angle a at Earth. The case k — +1 is shown 
in Fig. |. 

The objective is now to express the unknown side x'2 m 
terms of the known sides xi, X2 and the angle a. We first 
note that all formulae of Euclidean trigonometry have corre- 
sponding formulae in spherical and hyperbolical trigonom- 
etry. These can be expressed simultaneously for all three 
curvatures using E. In particular, we can generalize the half- 
angle formulae in this way and use them to show that 

v 2 fX2+Xi\ ■ 2 a 
— 2j sm 



(7) 



This is a more compact and symmetrical version of the gen- 
eralized cosine rule (see Section ^) . 

Using the same methods that were employed in the 
derivation of equation (^) we can relate the right-hand side 
of equation (R) to Z\ and z 2 : 
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p± 



where 
P+ = 

and 
P- 



1 

go 



H o1o v /(l + « 1 )(l + 2 ) 2 
(go - 1) (V 1 + 2q 2i + V 1 + 2 9o^2 - l) 



(8) 



+ y /{l + 2q z 1 )(l + 2q z 2 ) 
(^Jl + 2q z 2 - y/l + 2g 2i) . 



go 



(9) 
(10) 



Furthermore, setting xi = zi = it follows from equation 
(|) that 



aoE 



(¥) 



=== i (y / l + 2q z 2 -l 

+ Z2 2 \ 



Hoqo VI + Z2 
and for an observer at object 1 



(ii) 



(12) 



^igi yrr 

thus relating the left-hand side of equation (Q) to z' 2 . Solving 
P=\ J— r ( A /l + 2g^-l) (13) 

for 22 yields (positive solution) 
2P 2 



«2 = 



(«i " 2P 2 ) 



L + 'yi - ^/ >J + \/^ + i-2gi 



(14) 



where 



1 1+Zl 1 / 2 . 2 a 2 2«\ 

= T ^ ^ 7. P+ sm 7T ~\~ P— COS — (15) 

4 1 + 2 2 1 + 2g 2i V + 2 2 / ^ ; 

Note that z' 2 does not depend on Ho. 

2.1 A/0 

For completeness we now briefly consider the case A / 0. 
Equation (Q) is of course still valid. However, in relating this 
equation to the observables zi, z 2 and z' 2 we have made use of 
the fundamental relationship between an object's comoving 
radial coordinate and its redshift, 



ciqHo 



/l + z 
[u x 3 + (i - n )x 2 ] ' 



da;, 



(16) 



which in turn is derived from the Friedmann equations. For 
A / these take on a different structure and instead of the 
above relationship we have 



l + z 



[n x 3 + (l-Oo-A )i 2 + Ao] 2 di-,(17) 



where Ao = Ac 2 /3Hq. Unfortunately this integ ral is non- 
analytical but Kayser, Helbig, & Schramm (1997) developed 
an efficient method (which also accommodates inhomogene- 
ity) to compute E(x) numerically. 

Note that other density contributions with unusual 
equations of state can be dealt with in the same way. All 
that needs to be done is to establish x( z ) or i equivalently, 
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Figure 2. z' 2 as a function of a (angle on the sky between receiver 
and emitter), z\ = z 2 . Thick lines are for qo = 0.5, thin lines for 
50 = 0.15. 




Figure 3. z' 2 as a function of z\ = z 2 . Thick lines are for go = 0.5, 
thin lines for 50 = 0.15. 



£ X(^)]. In principle, one could then proceed to use equation 
([?) to find z 2 and r' 2 . However, when x( z ) i s n °t available 
analytically it is probably more practical to use the gen- 
eralized cosine rule (see Section ||, equation ^) instead of 
equation ([?]). In any case, one is still left with the problem 
of inverting £[^(2)] in order to find z' 2 . However, since we 
can compute 



(18) 



it should be possible to employ an efficient root finding al- 
gorithm for this task. 



3 RESULTS 

As cosmologists are often used to thinking in terms of red- 
shift rather than distance, we show the result of the above 
calculations in Figs. in terms of z' 2 , the redshift of object 
2 as seen by object 1. 

Figs, g and ^ explore the special case 21 =22. In Fig. ^ 
we show z' 2 (a) with z\ = z 2 fixed at various values. The thick 
lines show the case of a flat universe (f2o = 2go = 1) and the 
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thin lines show the case of an open universe with Qo = 0.3. 
Note that for large values of zi = z 2 there is some a™ such 
that z' 2 — > oo for a — > a M . This is the particle horizon 
of object 1 at the epoch corresponding to Z\. At that time, 
light emitted from objects separated from object 1 by angles 
> cioo has not had time to reach object 1 since the Big Bang 
(ignoring inflation) . In Fig. ^ we fix a at various angles and 
show z' 2 as a function of z\ = zi. The case a = 180° is 
often incorrectly used in undergraduate physics textbooks 
(e.g. [Halliday, Resnick, fc Walker 1993) , p. 1128-1129) as 
an example of how to add velocities in Special Relativity, a 
method which will give a wrong result for z 2 . 

In Figs. H and ^ we consider the special case a — and 
plot z 2 as a function of zi and z 2 respectively. Thick lines 
again represent a flat universe, thin lines an open universe. 
Whenever Z\ < z 2 (and a — 0), z' 2 is of course given by 



l + z' 2 = 



l + z 2 

1 + Zl 



(19) 



which is the only case where z 2 is independent of the cos- 
mological model, since in this case the time of emission of 
the photon received by object 1 is the same as the time 
corresponding to z 2 . Note that although XaOsfi, 22, a) = 
^2(22, zi, a) (cf. equation p3| ), the same does not hold for z' 2 
(cf. equations ^ and ^). Fig. |B| also provides the solution 
to an interesting thought experiment. What is the redshift, 
Zrcfl, of a photon emitted by ourselves (z 2 = 0) which was 
reflected back to us by a comoving mirror at z\l The answer 
is 



1 + Zrcfl = (1 + Zl)[l+ 22(21, 2 2 =0)]. 



(20) 



Finally we plot in Fig. g lines of constant z' 2 as a function 
of z 2 and a for a receiver at Z\ = 3. In this plot the flat and 
open cosmologies are represented by the solid and dashed 
lines respectively. 



4 RELATION TO OTHER SOLUTIONS 

We are aw are of two origi nal solution s to the probl em in the 



literature, Peacock (1999] ), p. 71, and Osmer (1981). Neither 



of them consider z' 2 . 
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Figure 6. z' 2 as a function of z 2 and a (emitter position). The 
receiver is at z\ = 3 and a = (marked by a diamond). The 
contours are lines of constant z' 2 . Solid lines are for go = 0.5, 
dashed lines for 50 = 0.15. 

Instead of equation (Q) , Peacock considers the general- 
ized cosine rule which directly gives E(x 2 )- Introducing the 
cosine equivalent of S(%), 



(21) 



(22) 



H( X ) = v/1 
we have 

S(X2) = H(xi) Sfos) + fc£(xi)£(xa) cos a 
which can be written as 

S 2 (X2) = E 2 ( Xl )H 2 ( X2 ) + E 2 ( X2 )H 2 ( Xl ) 
+ fcE 2 (xi)E 2 (x2) sin 2 a 
-2E(xi)£(x2)H(xi)H( X 2)cosa. 

When only r' 2 is needed the use of this equation seems more 
practical than our solution presented in Section ^. However, 
when using the gen eralized cosine rule the analogues of equa- 
tions (|l4| ) and (|l5|) are more complicated so that in cases 
where z' 2 is (also) of interest, e.g. when calculating the lu- 
minosity distance r' L2 = r' 2 (l + z' 2 ), the new solution is to be 
preferred. 



(23) 



On the cosmological distance and redshift between any two objects 5 




a (°) 



Figure 7. The contours are lines of constant , the frac- 



tional error made when using equation for z\ = 3 (marked 
by a diamond) and go = 0.15. The contour levels are indicated. 



Osmer's solution may be considered the most rigorous 
as it is based on a general resul t of differential geometry 
in maximally symmetric spaces. 



Weinberg (1972), p. 413 



showed how to transform to a coordinate system which has 
been 'quasitranslated'. Osmer uses this equation to trans- 
form from a coordinate system in which objects 1 and 2 
have particularly simple coordinates (see discussion in Sec- 
tion H) to one where the origin has been translated from 
Earth to object 1. The result is then given by 



E 2 (x2) sin 2 a 

+ [E( X a) 30a) cos a - E( X i) H( X2 )] 2 . (24) 



This equation of course reduces to equation ( J23| ) and thus 
the comments made there apply equally to Osmer's solution. 



5 VALIDITY OF APPROXIMATION 

Probably the m ost common approximation for the comoving 



dista nce is (e.g. Yoshii, Peterson, & Takahara 1993; Phillipp; 
1993) 



2 12 

a X2 



2/2 2v,2r /_m 2 

«oXA2 = [X(z)\ a 



+ 



a — (z)_ 



Az' 



where Az = zi — Z2 and z = (zi + Z2)/2 and 
a °Tz {z) = H®' 



(25) 



(26) 



In Fi g. u\ we plot the fractional error made when using equa- 
tion (Eq), 



apAx2 
aox'2 



aoXA2 - aoX2 
aox'2 



(27) 



against Z2 and a for object 1 at Z\ = 3 and for qo = 0.15. 
We see that the approximation gives both too large and too 
small distances (solid and dashed contours) depending on 
the position in the 2:2-0; plane. Note the 'ridges' along which 
the approximation incidentally gives the correct distance. 

The special case Az = deserves some further atten- 
tion as it corresponds to the well-known angular diameter 
problem: what is the length, L, of a rod at z which subtends 



an angle a as seen from Earth? Commonly, the answer is 
given as 

a cioE(x) 



La 



a r a 



1 + z 



(28) 



where rx is the angular diameter distance. However, strictly 
speaking La is the length of the line \ ~ const connect- 
ing the two ends of the rod at the epoch z, which is not 
the shortest distance between the ends. The length of the 
geodesic connecting the two ends is given by ax'2 which we 
can derive from equation (^). Setting \i = X2 = X an< i 
Z\ — Z2 = z we arrive at 



L = 



a-X2 

do 

~z 



2E" 



sinf 



(29) 



AL 
L 



The fractional error made when using equation (p8|), 
La j~ L , is of course the same as in (^B|) for Az — since 
the extra factor (1 + z) cancels out. However, comparing 
equations ( pj| ) and ( ]2s| ) it is particularly easy to see that 
approximations such as (^B|) actually contain two approxi- 
mations: (a) small angle and (b) neglect of curvature (for 
k = ±1). E.g. &k is independent of z in the flat case but is 
larger and varies with z for the open case. 
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